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Abstract. This article reviews recent progress of the analytical theory of quantum 
spin glasses (QSG). Exact results for infinite range and one loop renormalisation 
group calculations for finite range models of either insulating or metallic type are 
presented. We describe characteristics of fermionic spin glass transitions and of 
fermionic correlations which are affected by these transitions and by spin glass order. 
Connections between tricritical thermal- and T — QSG transitions are described. 
A general phase diagram with tricritical QSG transitions caused either by random 
chemical potential or by elastic electron scattering, and implying discontinuous 
T = O-transitions in weak and in strong filling regimes, is also derived. 

1 Introduction 

Fermionic quantum spin glasses form a part of the overlap regime of interact- 
ing disordered fermion systems and spin glasses. Understanding such systems 
with frustrated or with general random interactions between fermionic spins 
touches many currently active research fields. The wellknown transverse field 
Ising spin glass can also be represented as a pseudofermionic quantum spin 
glass with imaginary chemical potential and thus be treated analogously to 
the ones discussed here which involve real fermion spins. One of the attrac- 
tive features of fermionic spin glasses is the connection of spin and charge 
degrees of freedom. Even in a strongly localized Ising limit (with all spin 
and charge variables commuting) quantum statistics remains indispensable, 
since it governs the relative occupation of magnetic and nonmagnetic states. 
This provides a link between spin- and charge correlations, exerting a major 
influence on magnetic phase diagrams for example. Nonanalytic interplay is 
observed in the simplest fermionic spin glass, the SK-model extended to a 
four-state-per-site fermionic space (ISGf). This phenomenon occurs at a 
thermal tricritical point, separating continuous spin glass transitions around 
half-filling from discontinuous regimes at either strong or weak filling. This 
structure of the phase diagram is compared with a similar one of classical spin 
1 models including the clean BEG-model. Furthermore, model definitions on 
fermionic space generate a large set of correlation functions which display 
ubiquitous quantum dynamics of these models on the fermionic level, ie be- 
yond the one seen in spin- and charge-correlations and caused by electronic 
transport. The fermionic Green's function and its density of states content 
is one of those quantities which probe the commutation properties of single 
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fermion operators with the hamiltonian. The ISGf shows such quantum- 
dynamics. Moreover, correlations of the statistical density of states fluctua- 
tions turn out to be related to the spin glass order parameter, which intro- 
duces Parisi replica symmetry breaking into the charge sector. 



2 The Grand Canonical Fermionic Ising Spin Glass 

In contrast to standard spin glass models like SK- or transverse field mod- 
els, the models of interacting true fermion spins are naturally described 
in the grand canonical ensemble. Statistically fluctuating spin interactions 
leave only the choice between nonrandom local fermion concentration or 
nonrandom chemical potential. In any case however an effective spin dilu- 
tion due to thermal redistribution among magnetic and nonmagnetic states 
is controlled by either /i or v. This form of thermal spin dilution com- 
prises at a special point the case of nonanalytically communicating charge- 
and spin-fluctuations. The simplest model displaying this type of behav- 
ior is the Ising spin lass on a fermionic space (ISGf) with four states per 
site, defined by the hamiltonian H = —\ J2i^j Jij^i^j — ^J2i n i wrtn sphis 

a i = ^ r ta <T ap^ r i,l3 i particle number operator rn — ^j^La and gaussian 
distributed exchange integrals Jij [1]. The fermionic field operators obey 
the usual commutator relations {^q, ^73} = and {^,^73} = 8ij6 a fi- 
Expressing the partition function with the help of Grassman integrals the 
disorder averaging of the free energy is performed by means of the replica 
trick f3[F] av = [logZ] av = lim^oK 1 - [ n Z=i Za U)- As a guide to the 
global phase diagram we study an exactly solvable infinite range version of 
the model but the formulae obtained in this subsection may equally well be 
considered as saddle point approximation for an interaction with finite range. 



2.1 Phase Diagram 

Details of the calculation are given in [1], here we just state the result for the 
replica symmetry broken saddle point free energy 
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Notice the appearance of the Parisi variables q v and the additional q which 
lies at the heart of the following discussion. The Edwards - Anderson order 
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parameter qEA = lim^oo < Si(t)Si(0) > describes the freezing of spins in the 
spin glass phase and is given by the Plateau height q(l) of the order parameter 
function, whereas the replica diagonal q = [< a a a a >] av is related to the av- 
erage filling factor [u] av = ± J2ii n i\av by [is} av = l+tanh(/3/x)(l-<j). The last 
relation is exact even in the case of replica symmetry breaking. To obtain in- 
formation about the phase diagram a replica symmetric approximation is suf- 
ficient, though. The symmetric saddle-point solutions are q = J z smh c l^^ z ^ 

and q = jf cm gy ] with C^z) = cosh[/3H{z)} + cosh(/3 M ) cxp[-l/2/3 2 (g - 
q)] . Phase transitions are signalized by vanishing masses of the order param- 
eter propagators which in the saddle point formalism are given by second 
derivatives of the free energy. Hence a positive mass for q and a negative one 
for q guarantees stability. A similar system of coupled stability conditions 
was found for the BEG - model [2] and for a SK - model with crystal field [3] . 
Analyzing the stability limits one obtains curves of critical spin and charge 
fluctuations, respectively: 

Mcl (T)=Tcosh- 1 [(l/T-l)cxp[l/(2T)]] (2) 



Mc2 (T) = Tcosh- ^ 1 — exp[ i - Vi 2 _ == ]] (3) 

The two curves have a common tangent point at /i C 3 = 1/3 cosh -1 [2 exp(3/2)] = 
0.961056, T C 3 = 1/3. For smaller values of the chemical potential the param- 
agnetic solution becomes unstable to spin fluctuations first, at the common 
tangent point both types of fluctuations become critical simultaneously and 
tricritical behaviour results. 

2.2 Tricritical Point TCP: Exponents and Special Features 

An expansion of the saddle point equations around the tricritical point yields 
in leading order 

= gr g -r T 5T 2 + 66T5q- -5q 2 + 3q 2 , = 6q(Sq - ST - q) (4) 

where Sq = q — qTCP,gJ = A* — A*c3 + (C~ 1 J — Mc3)3<5T as nonordering field, 
and 5T = T - T c3 . The constants are given by r g = f^,r T = 2(1 - §C~ 2 ) 
with £ = tanh(n c z/T c z) ~ 0.9938. The average filling factor corresponding to 
/i c3 is evaluated as [^3]™ — 1.6625. From (4) we get for q=0 

Sq dls = 4(ST ± |<5T|^1 + T ^^g = 4(ST ± \ST\W) (5) 

Only the solution with the - sign corresponds to a minimum of the free energy 
in a region close to the line 5/j = —3(1/ C — ^cz)5T = — 0.1354<5T (tangent 
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to both ^ c i(T c3 ) and /x C 2(T , c3 )) g is of order ST 2 or smaller and usual critical 
behaviour results. However, if g is of order ST or larger the solution becomes 

to leading order Sq = ^ and thus displays a nonanalytical dependence 

on temperature and / or chemical potential. This type of crossover can also 
be seen from the scaling form of the free energy fdi S = \ST\ G(t§v) which 
allows for the identification of the specific heat exponent a = — 1 and the 
crossover exponent </> = 2. The crossover function Q{x) is regular for small 

x and has the asymptotic form Q(x) x ~ co (jj^) ? (<?oo+rcgular corrections). 
In the tricritical region the leading singularity in the free energy is given by 
fTCP.sing = "^( r 9-9)^- For <5/j = we have g ~ ST and can read off the 
tricritical specific heat exponent 013 = \ from above. 

Using q = Sq — ST one finds ordered solutions displaying the same crossover 
as discussed above, the transition from paramagnet to spin glass being con- 
tinous for positive ST and first order for negative temperature deviations. In 
the tricritical regime we find q = Sq = \\J—r g g which yields the tricritical 
order parameter exponent (3% — \ and suggests that q and Sq act as order 
parameters simultaneously. From the fluctuation Lagrangian (see section 3) 
one reads off mass squared proportional to ST and hence 73 = /?3 = as = |. 

3 Tricritical Landau Theory and the Parisi Solution of 
the Fermionic Ising Spin Glass 

We derived a fluctuation theory for the tricritical and finite range ISGf, for 
T ^ a Lagrangian of the same structure is obtained for models including 
transport mechanism by integrating out dynamical degrees of freedom 



l=\j d d X [ ^ £ + y: Q aa + \ e Q aa (-^ 2 + «)<?■ 



Here 4(^)^ = = (^)(°) - f and u (°) = denote the bare coeffi- 



cients at tricriticality. One fourth order term relevant for replica symmetry 
breaking is kept. Replicas under or Tr' are distinct. The Q oa Q hb- coupling 
is renormalization group generated as in the metallic quantum spin glass, its 
effects will be discussed in a subsequent section. 
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(d), starting at the TCP, 



tively. 



Inserting the Parisi Ansatz for replica symmetrie breaking in the Lagrangian 
(6) the expansion of the replica symmetry broken saddle point free energy 
around the tricritical point reads 



Here variables and constants are defined as in (3), r = (5T/T c3 . In contrast 
to crystal-field split spin glasses [3] the quartic coefficient yi of our free 
energy, Eq.(7), is nonzero and one obtains the Parisi solution q(x) = ^x for 
< x < x\ and q(x) = q(l) for x\ < x < 1. The plateau height is found to 
satisfy q(l) = Sq + 0(Sq 2 ). Consequently, plateau and breakpoint scale like 
\/b~l + 0(t) at the TCP, while linear r-dependence is reserved to T c > T c3 . 
Adapting the notation of [4] we express our result for the irreversible response 
q(l) - Jo q(x) ~ |t|^ in terms of the exponent 0A3 — 1 for T — > T C 3 and 
(3a = 2 for T — > T c > T c3 . For the Almeida-Thouless line at tricriticality we 



T c > T c3 . These values do not satisfy the scaling relation 83 — j^— with 

(3a 3 = 1 + (73 — a 3 )/2. Along the lines described in [4], this problem of 
mean field exponents will be resolved below by the renormalization group 
analysis of the coupling 3/4 of the finite-range and finite-dimensional ISG / . 

4 Replica Symmetry Breaking for Fermions 

The fcrmionic Ising spin glass allows for an exact evaluation of the bare 
fermion propagator G?y j(7 =<< - — 1 - » at least in the disordered phase; 





T c3 -T AT (H) 
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here H denotes the usual effective field and the double average refers to the 
replica-local and the Parisi block decoupling fields. The result can be written 
in the form 

r u\ - -i [J!! sr (2-a 2 M(^) 2 ^((i-a 2 )/3m.) , 8 , 



( e n - ip a + i\(3J 2 q ^ e» - ip a + i\j3 J 2 'q ^ 2 
( 6r}C{ Jy/2qJ ))£XP[ ^TqJ ] 6ih 



where fi a = fi + oH includes a magnetic field H and e„ = (2n + Vj-KksT/Ti. 
One easily extracts the disorder averaged electronic density of states [5] 

< P (e) >= -J— e -^ ch(f)ri+chWe + ri) 
^qJ ch{(3^)+ch{pH)e^ 2J2£ i 

Below the freezing temperature 0(<7 2 )-corrections occur in < p >. Statistical 
fluctuations Sp = p— < p > of the density of states are for example observable 
in < 5p®(e)5p®,(e') >. Taking Parisi symmetry breaking into account, this 
correlation becomes a function of the Parisi parameter x. Picking the one 
with = 0' and e = e' we obtain the (zero field) result 

< Sp a (e)Sp a (e) > (x) = A7tJ2 ~5 (1 - qc + q c e~?k cosh(^-^)) 2 

(1 + q c - ) 2 e-^q 2 (x) + 0(q 3 (x)). (10) 

This result shows that the fermionic density fluctuations reflect the irre- 
versible magnetic response introduced by the spin glass order. It is valid 
for all fillings, one only has to insert the appropriate Parisi function and the 
filling -dependent value of the spin autocorrelation q c at the critical point; this 
implies that the scaling of the density correlator with T — T c is different in 
the second order regime (quadratic) and at the tricritical point (linear), while 
for the discontinuous regime a one-step RSB is expected. It is also clear that 
for models with transport mechanism the calculation of conductance fluctu- 
ations is of great interest. So far we have put aside the question of replica 
symmetry breaking of fermion propagators: by this we mean the possibility 
of a nonvanishing propagator between different replicas, ie < ip ip b ^ a >. In 
the light of the Mezard-Parisi instability of the random field Ising model we 
feel that this preferably would occur as a fluctuation effect, if at all. 



5 Related Fermionic Spin Models 
5.1 The Fermionic Ising Chain 

Similarities between the phase diagrams of the clean BEG-model and of 
the fermionic Ising spin glass can be taken as indicative for the fact that 
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rather spin dilution than disorder is the source of the tricritical crossover 
from continuous to discontinuous phase transitions. While it is complicated 
to solve ID fermionic Ising spin glasses exactly, the clean fermionic Ising 
chain offers some simple exact solutions. Here we shall provide insight into 
the role of the chemical potential and moreover generalize known results into 
the complex /^-plane. Lee and Yang [7] derived the distribution of zeroes of 
the partition function of finite and infinite Ising chains within the complex 
magnetic field plane. Stimulated by the representation of conventional Ising 
chains by fermionic ones with special imaginary chemical potential, one may 
wish to extend the Yang Lee analysis to a fourdimcnsional space of complex 
(p, v). The transfer matrix Tf of the fermionic Ising chain reads 



T f = e fi » 



1 e -Pn e hP(h-») e -iP(iJ-+h) 

e ±P(p+h) e \p{h-y) e 0(J+h) e -0J 
\ e ±f)(p-h) g-i^+fe) e -f)J e f3(J-h) J 



(11) 



The transfer matrices Tf and Ts of the standard S — ±l-chain and their 
eigenvalues do not map onto each other at fj, = i^T, while the partition 

functions obey Z { f N) = TrT^ = (2i) N Z { S N) (fj, = i\T) for any number N of 
sites. The largest eigenvalue determines the free energy of the infinite chain, 
while the second largest is required in addition to determine the correlation 
length. The eigenvalues for h = are found as 

A± = e 0li [ch{(3n) + ch(j3J) ± y/{ch(fin) + ch(f3J)) 2 + 4c/i(/?/i)(l - ch((3J))} 
A = , Ai = 2e f ^sh((3J). (12) 

The correlation length is given by £ = l/ln(-^). In the T — > 0-limit a 
transition arises at fj, = J and due to the properties 

Thus £ diverges only for /x < J, since the energy for adding a fermion is larger 
than the gain from a magnetic bond if \i > J. Hence 

J exp((3(2J - fj,)) ,0<fi<J 
^\ exp(pJ/2) , M = J, ^ V 

while £ <~ T/(fi — J) for /i > J. The filling factor shows for T — > that 
the system is completely filled for fi > J (empty for fi < — J). Thus there is 
no physical T = O-transition of this simple system. The correlation length 
diverges in the T — > limit for all fillings v. The zero-field partition function 
shows that Yang-Lee zeroes approach fj, = ± J for T — > 0. This means that 
for (h = 0, T = 0) nonanalytic behaviour (as a function of the real chemical 
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potential) can only occur at the values fj, — ± J. It is instructive to consider 
N = 2 explicitly, which yields 

Zf = 2 = Ae 2 ^[{ch{l3fi) + ch{{3J)f + ch 2 {(ih){e 2f)J - 1) - ah(20J)] (15) 

This almost trivial case already shows zeroes at [i = ±(J + (^ln2 ± if )T), 
while allowing for finite complex magnetic field the first zero different from 
±J in the T — > 0-limit becomes possible with [i = ±(if + 2imn)T.More 
zeroes occur on the T — 0-axis as N is increased. For N — > oo a density 
function is expected in accordance with £ diverging for any [i. 



5.2 Mapping the twodimensional Ising model with imaginary 
magnetic field h = into fermionic space 

The complementary role of complex magnetic field and complex chemical 
potential can nicely be seen by recalling the exact solution for the 2d Ising 
model h = i^nT. This value corresponds to /j, = i^irT, which maps the 
fermionic Ising model onto the one above. Thus the exact solution of the 2d 
fermionic Ising model with \i = h = ij^T is known. Moreover this special 
model maps onto an interaction model of spinlcss fcrmions with a special 
species obeying bare Bose statistics but carrying along the minus signs of 
fermion interactions. The hamiltonian of this model can be written 

ij i i 

with fi = h = ijT. This reduces to zero the imaginary field of one fermionic 
species, while the other field equals the distance between Bose- and Fermi- 
Matsubara energies. Setting c~a^ and d — ai the hamiltonian reads 

H = - Jij \c\cic\cj - c\cid]dj - dldic]cj + d\did]dj] - d\di. (17) 

ij i 

The new imaginary chemical potential —inT of d-fermions renders their Mat- 
subara energies and single particle statistics bosonic. Thus anticommuting d- 
bosons interact with c-fermions. All vertices retain the fermionic minus-signs. 
In the fermionic path integral the c- and d-particles are described by the 
Grassmann fields Vc(t) and by Vd( T ) respectively. The i7rTn,i-term can be 
absorbed by the phase transformation exp(iirT)ipd(T) = ^(t). The new an- 
ticommuting fields ipd obey unusual bosonic periodicity ^d(r) = ipd(T + (3) 
on the imaginary time axis, while iP c (t) — —i\) c (j + 0) remains fermionic. 
The bosonic feature of the anticommuting fields Vd's is of course seen in 
the bare statistics. Perturbatively the above exact conclusion maps the 2D 
Ising model with h = inT/2 onto a coupled Fermi-Bose system (c-d) with 
additional (— l)-factors for each d-loop. 
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6 Metallic Quantum Spin Glass 

So far we have discussed the magnetic phase diagram of fermionic lattice gases 
having in mind that the tricritical phenomenon and universal quantities are 
not changed by coherent hopping of the electrons in flat bands. On the other 
hand it is well known [8], [9] that both a random chemical potential and 
an increasing band width suppress the transition temperature continuously 
down to zero thereby producing a quantum phase transition (QPT). The field 
theory of such a QPT 

Fig. 2. Phase 
diagram of fermionic spin glasses 
as function of chemical potential, 
hopping, and temperature. Above 
the marked area the system is in 
the paramagnetic, below in the 
spin glass phase. The shaded area 
corresponds to continuous transi- 
tion, the black area to first order 
transitions. Tricritical points are 
located at the boundary of the two 
regions. 



shows surprising similarities with the thermal tricritical theory, above all the 
simultaneous occurence of critical charge and spin fluctuations [9]. However, 
electron hopping not only lowers the T c of 2nd order transitions but also that 
of tricritical and 1st order transitions. For a metallic model with Gaussian 
random hopping and bandwidth 2Eq we found a quantum TCP at 

£ F = (l-(5/8)5)i£ , J c = 3^ [l-^/^o] _i /32, (18) 

Corrections to this Q - static approximation can be calculated by generalizing 
method [10]. The general phase diagram is sketched in figure 2. 

7 Renormalisation Group Analysis 
7.1 Tricritical Ising Spin Glass 

We performed a 1-loop RG calculation for tricritical fluctuations. At each 
RG step the mass of charge fluctuations 8Q aa was shifted away. Introducing 
the anomalous dimensions fj and r\ for diagonal and offdiagonal fluctuations 
one finds at one loop level the following RG relations (e = 8 — d) 

^ = - + 16kik 2 + &n\)r - k\, ^ = 2(1 - n\)u - + 4kik 2 , 
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ani e n an 2 e 22 aK 3 / e 2\ 

— = -Ki + 9k?, — = (- + 6k 2 - k x + 16kik 2 )k 2 , -^7- = + 9k 2 )k 3 , 

Above d — 8 the RG flows towards the Gaussian fixed point with mean field 
exponents, for d < 8 there is no perturbatively accessible fixed point for real 
k's. However, a preliminary analysis of the resulting strong coupling problem 
shows that there exists a solution with positive fj in contrast to the negative 
anomalous dimension typical of cubic field theories with imaginary coupling. 
The RG for the DIC x/a showed that its long-distance behaviour is dominated 
by a K 4 -contribution (like in [4] but) for d^ = 8 < d < 10. This leads to 
the modified MF exponent 8 3 — 8/(d— 4), which satisfies the scaling relation 
#3 = 2//3a 3 in d^ = 8 and reduces to the MF-result in 10 dimensions. The 
dimensional shift by 2 in comparison with [4] is due to coupling t. 

7.2 Metallic Quantum Spin Glass 

In the theory of quantum phase transitions time-dependent fluctuations are 
treated on an equal footing with spatial fluctuations [9]. While for the La- 
grangian (7) it was sufficient to consider only the lu = - component of the 
Q - fields, in the quantum case all low energy fields must be kept as they 
are coupled via the quantum mechanical interaction u J dr(Q aa ) 2 . However, 
the value z = 4 of the dynamical critical exponent renders the u - coupling 
dangerously irrelevant and allows for a perturbative mapping of the critical 
theory on a classical problem, the Pseudo - Yang - Lee edge singularity. This 
field theory has only one cubic coupling which corresponds to K\ in eq(7). 
The comparison of the metallic quantum case with the thermal tricritical 
theory allows one to understand the nature of the strong coupling RG fixed 
point: whereas the spin fluctuations in the thermal second order regime are 
governed by a perturbatively accessible fixed point, the TCP and the quan- 
tum case are charcterized by the combination of charge and spin fluctuations 
and a corresponding strong coupling problem. 
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